Fingerprints of the Magnetic Polaron in Nonequilibrium Electron Transport through 
a Quantum Wire Coupled to a Ferromagnetic Spin Chain 
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We study nonequilibrium quantum transport through a mesoscopic wire coupled via local ex- 
change to a ferromagnetic spin chain. Using the Keldysh formalism in the self-consistent Born 
approximation, we identify fingerprints of the magnetic polaron state formed by hybridization of 
electronic and magnon states. Because of its low decoherence rate, we find coherent transport sig- 
nals. Both elastic and inelastic peaks of the differential conductance are discussed as a function of 
external magnetic fields, the polarization of the leads and the electronic level spacing of the wire. 

PACS numbers: 73.23.-b, 75.10.Jm, 75.75.+a, 72.25.-b 



c3 



-a 
o 

(N 
> 

m 
m 
O 

o 

-i— > 

c3 



i 

C 

O 

o 



X 



Introduction. — In recent years, the field of Spintronics 
has attracted increasing interest 0, 0- A considerable 
amount of theoretical and experimental attention has 
been focused on transport phenomena, especially spin- 
dependent charge currents in low-dimensional structures 
made of magnetic materials 00,0], but also transport of 
magnetization through insulating spin chains and quan- 
tum dots 00. 

In this Letter, we study the interplay between nonequi- 
librium electron transport and magnetic degrees of free- 
dom in a one-dimensional system. The model under con- 
sideration is a finite quantum wire which is coupled via lo- 
cal exchange to a one-dimensional ferromagnetic Heisen- 
berg spin chain and via tunneling to two large electronic 
reservoirs. Examples of one-dimensional systems which 
exhibit ferromagnetic coupling of localized spins are so- 
called "sandwich clusters" formed from vanadium and 
benzene . Usually, one would expect that emis- 
sion of magnons in the spin chain will lead to a relax- 
ation and dephasing of the electron spins antiparallel to 
the spin direction of the spin chain, leading to incoher- 
ent transport for this spin direction. However, it was 
shown in several works on ferromagnetic semiconductors 
[ill Hi! Hii . Q] that a single electron with antiparallel 
spin direction to the localized spins can hybridize with 
one-magnon states to form the so-called magnetic polaron 
states. These states form a band which is separated from 
the band of scattering states, and therefore have a low 
decoherence rate. The aim of this Letter is to find finger- 
prints of these states in coherent transport signals at low 
temperatures by studying the differential conductance as 
a function of bias voltage. We note that one-electron 
scattering in finite quantum spin chains has been studied 
in Ref. 4] at low temperatures with the result of an in- 
teresting resonance structure as a function of the Fermi 
level. However, this work was restricted to linear trans- 
port, so that only states near the Fermi level contributed 
to transport. Therefore, the influence of the magnetic po- 
laron states (lying outside the band of scattering states) 
was not probed there. 

We calculate the differential conductance G = for 



large (N — 1000 sites) and small (N — 12 sites) sys- 
tems which differ in the electronic level spacing. We find 
peak structures which are due to elastic and inelastic 
transport processes. The applied magnetic field, the spin 
polarizations of the leads and the bias voltage affect the 
energies and decay rates of the electronic states of the 
system. One can thus control the position and height of 
the peaks in the differential conductance and identify the 
processes which contribute to the current. 

Hamiltonian. — We employ a tight-binding model for 
the quantum wire and a Heisenberg model with ferro- 
magnetic coupling J > for the spin chain (see Fig. EJ. 
They have the same lattice constant a and Zeeman split- 
tings h e and Ez, respectively: 



* E ( c U+i ff + h.c.) + -h e *4 



-ffspin — — J E S;Sj + i + Ez E Sf 



with % = 1, a = ± (we will frequently write a =\ or 
a =1 instead). We consider low electron densities in the 
wire and therefore neglect the Coulomb interaction. 

We follow [l2| and use the Holstein-Primakoff transfor- 
mation (HPT) 15] to replace the spin operators in the 
chain by boson operators: Sf » V2Sb\, 5r w y/2Sbi, 
Sf = b\b i — S. This approximation is valid if the spin 
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FIG. 1: The model under consideration. The conduction 
electrons can hop between neighboring sites in the wire, and 
the localized spins are coupled ferromagnetically. There is 
a local coupling between the spin of the conduction electrons 
and the localized spin at each site. The wire is coupled to two 
leads with the chemical potential ^il and fin, respectively. 
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chain is near its ferromagnetic ground state where (S 1 ?) = 
—S. The Zeeman energy Ez must be sufficiently large 
to ensure that this is the case. Using periodic bound- 
ary conditions, we can write H wilc = Y,ka £ ka c k<j c kcr' 
ff sp i„ = E + J2 k ^AK, where E = -NJS 2 - NSE Z 
is the ground state energy of the spin chain, b k and b k 
are creation and annihilation operators for magnons, and 
Ska = —2tcos(ka) + \oh e , uj k = 2JS(1 — cos(fca)) + E z 
are one-electron and one-magnon energies, respectively. 
If Ez ^> JS, we can assume that the magnon energies 
are independent of the wave number: lo^ lo = Ez- 

The interaction V = g^] t s^S^ between electron spins 
Si = \ J2aa> c L( <t )<j(j' c kt' anc ^ localized spins Si is trans- 
formed using the HPT to V = V w + V {2) + A, where 

kq 

corresponds to spin flips of a conduction electron which 
involve the emission or absorption of a magnon, 
m J2k qq 'a ah \+q> h q c \-q'o c ko implicates electron-magnon 
scattering, and A is a spin-dependent energy shift which 
can be combined with e k(T to form a new one-electron 
energy Ska = — 2t cos(fca) + \a{h e — gS). 

The wire is coupled to the leads by Hj- = 
J2 a ika (t>i la c ktT +H.c)). Here, a € {L,R} labels the 
lead and I the electronic states with energies e a i a . The 
leads are assumed to be noninteracting and to have a 
constant, but possibly spin-dependent, density of states. 
This is reflected in the energy-independent coupling func- 
tion = T%(E) = 27r£ ; t%t%*5(E - e alc ). 

We assume that the occupation of magnon states is 
equal to the equilibrium value n(u) = (exp(/3w) — 1)~ , 
i.e., that the coupling to an external spin bath causes 
magnon relaxation on a time scale tm which is smaller 
than the average time between two electron transmis- 
sions through the wire but larger than the time needed 
to establish a coherent electron-magnon state flif . 

Method. — The nonequilibrium Green function method 
proposed by Keldysh [16| is used to calculate the current 
through the wire which can be expressed in terms of the 
Green functions [17], We divide the electron self-energy 
into two parts, E = Et + Em, where Et is due to H^. Its 
retarded/advanced and Keldysh components are E^ = 
Tfr CT and = i £ Q T%(2f a (E) - 1), respectively, 

where r CT = ^ Q V£ and f a (E) is the Fermi function for 
lead a, which has the chemical potential fi a . 

The self-energy contribution Em, which is due to 
the electron-magnon interaction V, is calculated in self- 
consistent Born approximation (SCBA), which corre- 
sponds to the consideration of diagrams of the order 
0{g 2 ). We evaluate Em using the free magnon Green 
functions which do not depend on the magnon wave 
number due to the approximation LOk w lo. Therefore, 
Em is independent of the electron wave number. With 



M a {E) = J2 k \G ka (E)\ 2 , its imaginary part is given by 
ImE MCT (£ + ™) = -f^^j^rj-^), (1) 

a 

where /+ = / a , /~ = 1 — f a . The real part is obtained 
from the Kramers-Kronig relations. Terms which are pro- 
portional to n(uj) have been omitted because n(uj) <C 1 
for the parameters chosen in the next section. Dyson's 
equation and are solved self-consistently using an it- 
erative procedure [lfSj . 

Richmond 0] obtained the exact "["-electron self- 
energy for a single conduction electron in equilibrium by 
considering a larger set of diagrams involving an arbitrary 
number of electron-magnon scattering vertices between 
the emission and absorption vertices. However, these di- 
agrams would give rise to a violation of charge conser- 
vation in the nonequilibrium situation which is discussed 
here (i.e. the sum of the currents from the leads would 
be nonzero: Jl + 2r ^ 0). We therefore restrict our cal- 
culation to the charge-conserving SCBA and disregard 
diagrams of higher order than 0(g 2 ) [2fJ. 

Electrons which tunnel into the system from lead a 
with spin a can either tunnel to lead a' with unchanged 
spin a (elastic current) or flip their spin by emitting or 
absorbing a magnon and leave the system with spin —a 
(inelastic current). The two current contributions for 
lead a are 

T * = lj dE EE T ««v(£) (/«(^) - few > 

a' (T 

C C V) = lfdEj2 (.TZ}(E,Lv)f a (E) [1 - f a ,(E to)] 

a' ^ 

+ T^{E,u)[l-f a {E-u)]f a ,(E)y 

where we have again neglected terms cx n{u>). The trans- 
mission coefficients are given by 

T*,„(E) = T^M a (E), (2) 

tZ\e,lo) = 0r?rfM T (£)M;(i;-u,). (3) 

We consider both nonmagnetic and ferromagnetic leads 
with polarization P a = (1^ - Tf)/(T^ + T%) ^ 0. 

Results. — The differential conductance for a large sys- 
tem with N = 1000 sites where the level spacing is 
smaller than the other relevant energy scales is shown 
in Fig. |2 Since the electron density in the wire should 
be sufficiently low to justify the neglect of the Coulomb 
interaction, we examine the voltage regime where only 
a small fraction of the t - and |-electron states is par- 
tially occupied. We fix /^r below the conduction band 
at /xr = — 2.5 and vary (all energies are in units of t 
which we set to t = 1). 

For low bias voltages (/iL ^ —1.85), elastic transport 
processes (corresponding to resonances of |G^ (7 (^)| 2 ) 
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FIG. 2: Differential conductance for a system with N = 1000 
sites and S = 1/2. The parameters are t = 1, g = 0.5, 

E z = h e = o.i, k B T = 5xht 4 , £ CT r^' R = i0" 2 , p l = +0.7, 

Pr, = —0.7. The chemical potential //r of the right lead is 
fixed at —2.5. Inset: Conductance for unpolarized leads (solid 
line) and ferromagnetic leads with parallel polarization and 
P = +0.7 (dotted line), P = -0.7 (dashed line). 



FIG. 3: Conductance for N = 12, S = 1/2. The parameters 
are * = 1, g = 0.1, E z = h e = 0.04, k B T = 5 x 10~ 4 , 
I] CT r CT' R = 2 x 10 ~ 3 , MR = -2.5. Both nonmagnetic leads 
(solid line) and spin-polarized ferromagnetic leads with Pl = 
+0.7, Pa = —0.7 (dashed line) are considered. Inset: G for 
ferromagnetic leads with parallel polarization. 



dominate. One peak at fxi, ~ —2.12 for spin-up is 
due to the magnetic polaron band, and the other at 
/il ~ —1-89 corresponds to |-electron states. The shape 
of the peak structures reflects the density of states in a 
one-dimensional system. Inelastic processes superpose 
these structures. In the situation where Pl = +0.7 
and Pr = —0.7 (main plot of Fig. they dominate 
for higher voltages (fij, > —1.85) because this setup 
maximizes the product T^r^ which is proportional to 
the inelastic transmission coefficient The inset of 

Fig. [21 shows how the different peaks can be distin- 
guished by changing the spin polarization of the leads. 
Here, three conductance curves are shown for unpolarized 
leads and ferromagnetic leads with parallel polarizations 
(Pl = Pr = P). The weight of the peak structure which 
is related to elastic transport of electrons with spin a is 
greatest if the leads are a-polarized (a e {f , J,}), just as 
one would expect. While there is also an inelastic current 
contribution in these configurations, it does not cause a 
clear signal in the differential conductance. 

For a small system with TV = 12 sites, the discrete 
structure of the energy spectrum can be identified in the 
differential conductance if the level spacing is larger than 
the energy scales T and fcsT which determine the broad- 
ening of the conductance peaks. In the considered voltage 
regime, only the lowest electronic states (wave number 
k = 0, spin f or J.) are partially occupied and contribute 
to the current. The electronic spin- up states are split 
mainly by the q = magnon into two states with en- 
ergy £| and (corresponding to the magnetic polaron 
and the scattering state in the continuum case, respec- 
tively). Futhermore, the electronic spin down states are 
also renormalized to when the f- electron states have 



a finite occupation probability. 

Results for g > (antiferromagnetic local exchange 
coupling) are presented in Fig. [3J As in the situation 
discussed above, conductance peaks which are due to 
elastic transport processes coincide with resonances of 
the retarded Green functions. Here, the left and right 
large peak occur at £j and £|, the main resonances of 
the j- and |-electron Green function (with wave num- 
ber k = 0), respectively, and can be attributed to elastic 
transport through the wire. On the other hand, the in- 
elastic transmission coefficient (0 is proportional to both 
|Gj?|(P)| 2 and |Gj^(P — w)\ 2 . Therefore, inelastic trans- 
port processes contribute to the differential conductance 
at ^l = £| ( wnere some weight is added to the left large 
peak) and /Ul = +u>, the position of the smaller peak. 
The dependence of the relative peak heights on the lead 
polarizations is like in the large system discussed above. 

Actually, one could expect another peak in Fig. be- 
cause the t-electron Green function has a second reso- 
nance at an energy ei. However, it has in general quite 
a small weight because the magnitude of the imaginary 
part of the t-electron self-energy is rather large at et, 
leading to a strong decay of the corresponding state and 
a suppression of elastic transport. Moreover, is very 
close to the energy e j + uj where inelastic transport con- 
tributes to the current. Therefore, elastic transport of 
t electrons is visible only at £| (energy of the magnetic 
polaron), the small contribution at /il = £j is absorbed 
in the inelastic peak. 

It should be noted that not only the peak heights but 
also the positions depend on the polarizations of the 
leads. Fig. shows that the peak at ^l = ej. and the 
inelastic peak are shifted to the right for the configura- 
tion where the leads have antiparallel polarizations. The 
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FIG. 4: Conductance for t = 1, g = -0.1, k B T = 10" 3 , 
F^ R = 5.5 x 10~ 3 , fin = -2.5. The Zeeman splittings h e 
for conduction electrons and Ez = 2h e for localised spins are 
different for each curve. 



reason is that a J, electron (or hole) can interact with 
magnons only by nipping its spin and occupying an | 
state. Therefore, the |-electron self-energy (and thus 
the position £j of the main resonance of the Green func- 
tion) depends on the occupation probability of |-electron 

pa 

states. This probabilty is given by F^(E) « ^ Q j^f a (E) 
for a state with energy E and is affected by both the 
chemical potentials and the polarizations of the leads. 

Conductance curves for ferromagnetic local exchange 
coupling (g < 0) and different magnetic fields are shown 
in Fig.0] The Zeeman splitting Ez of the localized spins 
is chosen to be twice as large as the splitting h e of the 
conduction electrons. Therefore, not only the peak po- 
sitions but also the general structure of the conductance 
curve change if the magnetic field is varied. For small 
(h e « 0.02) and large (h e « 0.06) fields, the situation is 
comparable to Fig. [3] There are two large "elastic" peaks 
at /il = £f and /iL = £j. and a small "inelastic" peak at 
/iL = ej, +u>. These peaks move with different 'velocities' 
if the field is increased: The positions of the large peaks, 
i.e., of the main resonances of the retarded Green func- 
tions, change with the conduction electron Zeeman en- 
ergy ±h e /2, but the position of the inelastic peak changes 
like —h e /2 + Ez = 3/2h e because of our choice Ez = 2h e . 
One could expect that the inelastic and the right elastic 
peak overlap and form a single resonance for intermediate 
fields (h e w 0.045), but this is not the case. The corre- 
sponding conductance curve rather reveals two peaks of 
comparable height. These arise from two resonances of 
the "["-electron Green function which have approximately 
equal weight for this particular set of parameters. This 
means that the decoherence rates are equal for the states 
corresponding to the energies e-\ and e*, in contrast to 
the situation in Fig. Elastic transport of | electrons 
thus generates a double-peak structure in the differen- 
tial conductance which is superposed by a small inelastic 
transport contribution. 



We remark that without coupling to the spin chain 
(g = 0), there would be only elastic transport through the 
wire. The effects discussed here, i.e., inelastic transport, 
a dependence of peak positions on the lead polarization, 
and magnetic field-dependent decoherence rates of the 
states involved in transport, would not occur. 

Summary. — We presented a self-consistent diagram- 
matic approach within the Keldysh formalism to cal- 
culate the nonequilibrium current through a mesoscopic 
quantum wire coupled to a ferromagnetic spin chain. We 
proposed a way to detect the coherent superposition of 
electronic and magnon states, the so-called magnetic po- 
laron. It shows up as a high (i.e. coherent) signal in the 
differential conductance and can be tuned by external 
magnetic fields and the spin polarization in the leads. In 
this way we have shown that the interaction between elec- 
trons and magnons (which usually leads to unwanted re- 
laxation and dephasing of the electron spin) can be used 
for the creation of a phase-coherent quantum state. We 
expect that this work will stimulate further theoretical 
and experimental investigations of the magnetic polaron 
in the field of mesoscopic systems. 
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